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Abstract 

We study the Be meson rare decay in order to search for the Majorana neutrino signal. It is 
found that the the corresponding decay rate is sensitive to the Majorana neutrino mass and mixing 
angles. The signal of B^ — > l^l^M^ induced by the Majorana neutrino within the mass region 
m n < m n < mB may be observed at LHCb. 
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I. INTRODUCTION 



In the Standard Model (SM), the neutrinos are massless since there is no right-hand 
states. However, non-degenerated neutrino masses provide the most accepted explanation 
for neutrino experiments [lMSl • This is the favorite evidence for new physics beyond SM. The 
neutrino masses can be obtained by including the right-hand states, just like the treatment 
for all other fermions via Yukawa couplings with the Higgs doublet in SM. But unnaturally 
the hierarchy problem will become more serious due to the small mass of neutrinos. Other- 
wise as the right-hand neutrinos are SM gauge singlets, the Majorana mass term cannot be 
ruled out by the gauge invariance. In fact under the help of the Majorana mass term, it could 
naturally explain the smallness of the neutrino mass with the so-called see-saw mechanism 
[^l . The particular interest in this regard is the question about whether neutrinos are Dirac 
or Majorana particles. A crucial role to address this question is being played by several 
experiments looking for the possible existence of the lepton number violation processes. 

There have been several attempts to determine the Majorana nature of neutrinos by 
studying the lepton number violation processes. The experimental observation of such pro- 
cesses may be induced by the Majorana nature of neutrinos. The neutrinoless double beta 
decays (0vf3(3) in nuclei, (A, Z) — > (A, Z + 2) + 2e~, have been studied widely. By assuming 
that Oz//3/3 in nuclei are mediated by the exchange of light Majorana neutrinos, the higher 
precision of present experimental data has been able to set strong constrains on the effective 
mass, (m n ft) = J2 V U av Up v m Vl where a,(3 = e, fi, r. The upper bound on (m ee ) is 0.2 evficl]. 



A global fit 



ll| at 99% C.L. gives 1.1 x 10 3 eV < (m ee ) < 4.5 x 10 3 eV A (normal hierarchy) or 



1.2 x 10 3 eV < (m ee ) < 5.7 x 10 3 eV (inverted hierarchy). From the analysis of atmospheric 



and solar neutrino oscillation and the tritium beta decay endpoint experiment [12J, [13| , the 
limit is (m^) < 4.4eV^[14]. However, it has long been recognized that, even though the 
experiments are very sensitive, the extraction of the properties of the Majorana neutrinos 
from nuclear Qv{3{3 is a difficult task, because it is reliable only if the nuclear matrix elements 
for Qv/3/3 are calculated precisely. 

Another way to detect the Majorana nature is to study the lepton number violation 



processes pp — > l^l^ 1 + X at LHC 15l4l8l|. The lepton number violation processes in meson 



rare decays have also been investigated in refs. 



19h27|]. The aim of this work is to investigate 



the Bq — > UM. This signal process can be captured at high intensity experiments such as 
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LHCb 



281 ] and future super B factories. The AL = 2 processes B c — > I I M T can occur 



via Majorana neutrino exchange, and thus their experimental observation is helpful to test 
the Majorana nature of the neutrinos. 

The paper is organized as follows. The Lagrangian of Majorana neutrinos is introduced 
in Section [Til The formulas of the Be —> 11 M decays are obtained in section IIHI In Section 
IIVI we give the numerical results and discussions. Finally, a short summary is given. 



II. LAGRANGIAN RELATED TO MAJORANA NEUTRINO 



With the same gauge group SU{2) L <g> U(l)y in SM, the leptonic content in the simplest 
extension of the SM includes three generations of left-hand SU (2) L doublets and n right- 
hand singlets 

v 
I 

One can write the general gauge invariant Yukawa terms with Majoranan mass terms of 
right-hand neutrinos as 



/ 

V 



l R , Nr. 



(1) 



L Y = fiL L $l R + f u L L $N R + NrMrNr + h.c. 



(2) 



Therefore, the complete neutrino mass sector is composed of both Dirac mass which is 
produced via the Yukawa couplings with the Higgs doublet in the SM, and heavy Majorana 
mass term. 



-C 



M 



u L M D N R + N c rMrN r + h.c 



{"l,Nr) 



t 



R 1 

M D 

T 



Mr 

One can find the light Majorana neutrino mass are 

M v ~ -M d Mr x MI, 



N, 



+ h.c. 



(3) 



R 



(4) 



which is called Type I see-saw mechanism. There are other proposals to naturally generate 
Majorana mass for neutrinos called Type II or Type III. Generally the mass terms of neutri- 
nos with both Dirac and Majorana terms after gauge symmetry breaking can be expressed 
as 



- C = i7 L M D N R + ^N C R M R N R + l -v L M L v c L + h.c. 



(5) 



M 



To diagnose the mass matrix 

Ml M R 

one need introduce mixing matrix between the gauge and mass eigenstates 

/ 



T 



3x3 V3 Xn 



\ ^nx3 ^iixn 



The mixing matrix is unitary, TT^ = T^T = 1. The states are redefined as follows, 



(6) 



(7) 



ViL -> UijUjL + V ik Nt 



kRi 



n r -)• Xyi/g + f^aw 



(8) 
(9) 



As the M N > M v , one has 1/ ~ X ~ M D M^ and WM R + Y^M D ~ 0. 

In terms of the mass eigenstates, the gauge interaction Lagrangian of the charged currents 
now has the following form, 

„ / r 3 t 3+n \ 

£ = -4^7 E E utn&zfPrf + E E vs^Vw + ^- (10) 

where Pl = |(1 — 75), v m {m = 1, 2, 3) and N m r(mf — 4, • - • , 3 + n) are the mass eigenstates, 
t/^ m is the mixing matrix between the light flavor and light neutrinos, and V£ m / is the mixing 
matrix between the light flavor and heavy neutrinos. 



III. DECAY WIDTH FOR B% -> iflfM* 




The Feynman diagrams which contribute to Bq — > lflfM T are shown in Fig.([T]). The 
first diagram is the annihilation which has been studied widely with the input of meson decay 
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constants. The second one is the emission, which is considered not enough due to the non- 
perturbative long distance contributions in the hadronic matrix elements. The amplitude of 
B c decay to pseudo-scalar contributed by annihilation diagram is expressed as 

M% = -i2G 2 F V B V M f Bc f M (yCp ■ 1Ps ■ 1 P L v 2 ) £ S„(p a ) + (1 <+ 2), (11) 

V 

where Vb = Vcb, and Vm is the corresponding CKM element for final meson, p is the 
momentum of initial Be meson. Mb is the mass of the meson Be- P3 is the final meson 
momentum and Pxjpi is the momentum of lepton. The /# and ju are the decay constants. 
For a charged or neutral pseudo-scalar meson, the decay constant is defined as 

(OIAJP*) = if P p„ (12) 
V2(0\A,\P°) = if pPfl , (13) 

where the operator stands the axial-vector current operator. 
The full expression of propagator function S„ is 

c t \ V hu Vi 2U m u 

SM = (14) 



The V\ iV is the lepton mixing matrix element. With the same denotation, the amplitude for 
Bq decay to vector (or scalar) is 

My = i2G 2 F M v V B V v f B f v {vgp ■ 7 e 3 • jP L v h ) £ S v ( Pl ) + (1 2), (15) 

M% = -i2G F V B V M f Bc fM mi ~ m2 (vfp ■ IPs ■ -yP L v 2 ) ]T S„( Pl ) + (1 o 2). (16) 

In addition to the light neutrinos, the see-saw mechanism predicts very heavy neutrino also. 
The propagator functions for the two type neutrinos can be approximately expressed as 

S(pi)^{ ^ +P v f (17) 
{ Yh ^ L m N »M B , 

where we still use v to denote for the light neutrino, and N for the one heavier than tub- 

In the meantime, there is another possibility that the neutrino mass is between the Mb 
and the final meson mass which is denoted with n (m^ < m n < m B )- Neutrinos with such 



mass has been strongly constrained by direct search and cosmological observations 29N32 1 
and must be sterile. The process of Bf — > lfl 2 M T is dominated by the annihilation diagram 
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as the intermediated neutrino in s-channel can be on-shell. At this case the narrow width 
approximation 

1 7T 

lim — — — — — = — —8(q 2 — m 2 ), (18) 

r^o(g2_ m 2)2 + m 2 r 2 mT vy h \ ) 

can be applied. As the + X are the dominated decay channels, the total decay width of 
the neutrino can be expressed asfl9|. 

r„ = 2]T^ n | 2 (^) 5 xr r , (19) 

where the m T , T T are the mass and total width of the tau-lepton. With the total decay 
width and narrow width approximation, one can get the decay width of Bq — > lfl 2 n T 



ml I \ ml 



In general, the square of annihilation amplitudes can be written as 

\M a P \ 2 = 8(G 2 F V B V M fBcfM) 2 (pi-P2)F(p u p 2 ) + ( Pl ^p 2 ), (21) 



\M a s \ 2 = 8(G 2 F V B V M f Bc f M ) 2 ( mi Ms m2 ) (pi-P2)F( Pu p 2 ) + ( Pl ^p 2 ), (22) 
\M a v \ 2 = 4 (G 2 F V B V M f Bc f M f My pPh^ pfa ) + (pi „ ) (23) 

where the function F(p 1 ,p 2 ) and V(pi,p 2 ) are defined as, 



F{ P i,p 2 ) 



\Y, v V llv Vi 2U m v \*, m I ,<m 7r , 
\V hn V l2n \ 2 ^5((p - Pl ) 2 - ml), on shell, (24) 
\V hN Vi 2N \ 2 ^£, m N ^>m B , 



2 , 4(pi ■ P3 f 



V( P x,p 2 ) = 2pi • p 2 ^4pi -p 3 - M v -\ — 2 J + 8pi ■ p 3 p 2 ■ p 3 . (25) 

In addition to the annihilation diagrams, we also consider the emission diagrams. We 
apply the light-cone function of mesons to calculate the hadronic amplitude. The non- 
perturbative effect from long-distance interaction is factorized into the light-cone function, 
and the leptonic number violation effect is caused by the short-distance interaction which 
can be calculated perturbatively. The light-cone distribution amplitude is defined as {3^ 

r rl^h 



Mf a {k) = I j^e-^{0\q a (x)qp(0)\M), (26) 
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and parameterized with the twist wave functions as 
,b _ ifB 



{(^ + m B ) 7 5 B ( M )}, 



(27) 



W = 4UW-M 5 ( M«) ~ i^n^^P- + ia^jf ^ J )\. (28) 
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6 dk 



Tu 



35|, 



For the wave function of B c , we take the following form in the numerical calculations 34 



1 / X1TLb\ 2 

2 \ 



(29) 



4>b{x) = Nbx 2 {1 — x) 2 exp 

where ub is the shape parameter, and Nb is the normalization constant. The general 
expression of twist-2 wave function for pion is 



4>k{ x ) — — x) 



1 + J2 *nCl /2 (2x - l)\ , 
\ n=l / 



(30) 



where C n (x) is the Gegenbauer polynomial. In this work, the higher twist contributions are 
not considered. 

The amplitude of the contribution from emission Feynman diagrams can be written as 



M 1 



dxdy 2G 2 F V ub V cd 



{™>hi 3 ) _|_ ViinVi 2 n 



(31) 



where q = xps — ypz — Vi- x and y are the momentum fractions. The first term in Eq. (!3T]) 
is the contribution of light neutrinos while the second term is that of heavy neutrinos. The 
leptonic and hadronic tensors are 

(32) 



= (vr| (V(l - 7 5 )m) (J 7 v (1 - 7 5 )c) \B C ) 

= W« + p v rt - p ■ p*!T) MWM, 

Zl\ c 

where Nq is the color factor. One can get the emission amplitude as 

E 2G 2 F V ub V cd f B f 7r P ■ Pn<j>B<Pir 



(33) 



M 1 



x 



dxdy 



N c 
( m hh) 



+ E 



VhNVl 2 N 



+ (1«2). 



(34) 



q z + le — m N 

Comparing Eq. (|34p with Eq.f llip . one can notice that the contribution of heavy neutrino in 
emission diagrams is similar to that in annihilation diagrams but suppressed by the color 
factor Nc, while in some channel also suppressed by the CKM matrix element. Thus the 
emission diagrams can be considered only for light neutrinos. 



p 


7T (ud) 


K~(us) 


D-(cd) 


D-(cs) 


B-(ub) 


B+(cb) 


fp 


130.41 


156.1 


206 


257.5 


193 


322 


Mp 


139.57 


493.677 


1869.60 


1968.47 


5279.17 


6.277 



TABLE I: The input parameters for pseudo-scalars 



V 


p (ud) 


K*-(us) 


D*-(cd) 


Dries) 


fv 


216 


220 


240 


272 


my 


770 


891.66 


2010.22 


2112.3 



TABLE II: The input parameters for vector mesons. 
IV. NUMERICAL RESULTS AND DISCUSSION 



To compare the contribution from the annihilation diagram with that from emission one, 
we study them seperately. The numerical values of the mesons' decay constants and mass 
are listed in Table. ( ]T][TT|) . For the light neutrinos(m^ <C m n ), we take the effective 



36 



mass of the light neutrinos (m^ 2 ) = \J2vVhvVi 2 v m v\ as input. The contribution from the 
annihilation diagrams is found to be 



Br(Bf lflfn T ) ~ 1.4 x 10" 



3.3 



leV 



One can find the results for other channels in Tabel. (IHIj) with (mi^) = leV. 



p 


7T (ud) 


K-(us) 


D-(cd) 


D~(C3) 


Br 


1.4 x 10~ 33 


9.6 x 10~ 35 


6.5 x 10~ 35 


1.8 x 10~ 33 


V 


p~(ud) 


K*~(us) 


D*-(cd) 


Dries) 


Br 


6.0 x 10~ 32 


2.3 x 10~ 33 


3.0 x 10~ 34 


6.4 x 10~ 33 



(35) 



TABLE III: The branch ratio of B+ — > ifl^M with (mi^) = leV from annilation. 



The results indicate that the brach ratio of Bf — > iflfP^ is about 10~ 33 x |Vp| 2 x 

{ ^"lev^ ') ' ^ )ue ^° ^ ne su PP ress i° n °f f ne CKM elements V us ,V c d, the branch ratio to K(K*) 
and D(D*) are about one order smaller than n(p) and D S (D*). As the possible relevant phase 
between V^and Vi 2U , the (mi lU i 2U ) may be very small which indicates invisible contributions. 
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In addition to the annihilation diagrams, the contribution from emission diagram is 
depend on the non-perturbative parameter ujb- For lob = l.OGeV, the branch ratio of 
Bf — > / ± / ± 7r =F contributed by emission diagram is 

Br(Bf Iflf^) = 6.75 x 1(T 22 {^^j ' • (36) 

In Fig.(T5]) the result as a function of ujb is shown. One can notice the emission diagram is 
dominated comparing to the annihilation diagrams for light Majorana neutrinos. However 
it is still below the experimental bounds. 

io- 19 

10 -20 

io- 21 

IO" 22 
IO" 23 

1 2 3 4 5 

w Bc GeV 

FIG. 2: The branch ratio of — > l^l^ir^ as a function of the shape parameter ujbc- 

Now we consider the effect of heavy neutrino(m7v > m^). The nature explaining to the 
smallness of neutrino masses in the see-saw mechanism need 3> Aew where Aew is 
the Electro- Weak scale. The mixing element Vin is also very small. One can image that 
V1NV2N / mN ~ 1/Agi/t- If ^ijv^2Jv/wjv = 10~ 16 GeV~ 1 is taken, one can get the annihilation 
result 

Br(B- -» HH^) K 3.3 x 10~" ( 'i™^"" )'- < 37 > 

The results related to other channels with ViatT^jv/^jv = 10 _16 Ge^ _1 are listed in Ta- 
ble. fllVI) . One can notice that the contribution of such massive neutrinos to B^ — > lflfM T 
can be omitted comparing with contribution of light neutrinos. As pointed in last section, 
the contribution of massive neutrinos in emission diagrams is suppressed and more less. 
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p 


7T (lid) 


K-(us) 


D-(cd) 


D-(cs) 


Br 


3.3 x 1(T 45 


2.1 x 1(T 46 


9.6 x 1(T 47 


2.6 x IO" 45 


V 


p~(ud) 


K*-(us) 


D*-(cd) 


D* s ~(cs) 


Br 


6.1 x 1(T 45 


3.4 x 10~ 45 


3.5 x 1(T 46 


6.6 x IO" 45 



TABLE IV: The branch ratio of -»• ifl^M^ with V 1N V2 N /m N = 10- 16 GeV -1 . Only the 
anihilation diagrams are considered, since the emission diagrams are suppressed. 





m v « Tn^ 




m n » tub 


Emission 


~ io- 21 


~ io- 21 


< IO" 46 


Annihilation 


~ 10~ 32 


< IO" 4 


~ IO" 45 



TABLE V: The comparing of the contribution from emission and annihilation diagrams. 

At last we discuss the sterile neutrinos with mass m v < m n < tub- We take m T = 
l.llGeV, T T = 2.3 x 10~ 12 GeV and m n = ms/2. The intermediated neutrino in annihilation 
diagram can be on shell which will enhance such process. The emission diagram contribution 
has no enhancement which can be omitted. We list the approximational contribution from 
emission and annihilation diagrams in Table. El With such inputs, the branch ratio as 
functions of the leptonic mixing matrix elements are shown in Fig.flS]). The numerical result 
of branch ratio for Bq — y lflfn T is 

Br(B+ -> Iflt^) ~ 4.39 x IO" 4 x Ili!^!!. ( 38 ) 

Hi \ Vln\ 

The J± n , ' 2 " 2 is a parameter relating with mixing between the leptons. The brach rations 

/ , i I 'in I 

as a function of the parameter are shown in Fig. (J3j) . One can notice that the rare decay 
7r T may be detected on LHCb at this case. It is expected to study the property 
of the majorana neutrinos indirectly through such processes. 

The production of B c at LHCb have been studied widely. The gluon-gluon fusion sub- 
process gg — > Be + X is the dominated production channel and much larger than the 
quark- ant iquark annihilation subprocess qq — » Be + X. The magnitude of the color-octet 
components may be estimated with the non-relativistic QCD and the contribution can be 
neglected comparing with color-singlet components. 
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As shown in Ref 



40], the corresponding cross-section is about 3nb at Tevatron(-y/i 



1.96TeV) and 50nb at LHCb(y / i = 14TeV). As the desired Luminosity of LHCb is about 
10/6 -1 , the expect event numbers are also shown in Fig. ([3]). The LHC has run at 7TeV 
with total Luminosity and the production cross-section of B c is about 22nb. The 

exclude region has been shown in Fig. (HJ) at 95% C.L.. 




FIG. 3: The branch ratio of Bq and the corresponding event number on LHCb as functions of 
^"p^jjl with rn n = m B /2. 



i 

o.i 
o.oi 
E 3 o.ooi 
io- 4 
io- 5 

10~ 6 




FIG. 4: The solid lines stand the low bound of the parameters that is excluded by Tevatron and 
LHCb through B^ — > lflfir T at 95% C.L.. And the dashed line is expected to be excluded by 
LHCb in future. 
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V. SUMMARY 



The neutrino experiments indicate very small neutrino mass. A most natural explaination 
to such small mass is see-saw mechanism. In this work, we study the AL = 2 semi-leptonic 
decays of Be meson mediated by Majorana neutrinos and investigate the contributions from 
annihilation and emission diagrams for different Majorana neutrino mass. The light-cone 
functions of the mesons are applied to calculate the hadronic matrix element in emission 
diagrams. It is found that the corresponding decay widths are sensitive to the Majorana 
neutrino mass and the mixing angles. For a sterile neutrino with mass m n < m n < tub, the 
leptonic number violating decay rates of B c can be enhanced by the annihilation diagrams 
and may be detectable at LHCb. 
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